Electronic and magnetic properties in strongly correlated heterostructures 
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We present a theoretical study of a model heterostructure for a Mott-insulator sandwiched between 
two band insulators, such as SrTiOs/LaTiOs. Particular emphasis is given on the interplay between 
magnetism and inhomogeneous charge distributions. Our mean-field analysis of the generalized 
Hubbard model displays numerous ordered phases in the ground-state phase diagram. In particular, 
we find a canted antiferromagnetic state near the interface when antiferromagnetic-ordering exists 
inside the Mott insulator. A checkerboard charge-ordering proposed previously is also stabilized for 
large long-range Coulomb interactions. Regarding its origin we also point out the importance of 
interlayer spin-mediated interactions. It is further shown that such a strong spin-charge coupling 
gives rise to pronounced magnetic/charge order phase transitions in external magnetic fields: a first- 
order metamagnetic transition and a reentrant charge-order transition with checkerboard pattern. 
The mechanisms stabilizing these intriguing phases are explored through a detailed analysis of the 
physical quantities with special focus on the spin-charge interplay. 
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I. INTRODUCTION 

Correlated heterostructures composed of transition- 
metal oxides have attracted much interest both exper- 
imentally and theoretically as a new arena for studies 
of strongly correlated electron systems. Especially, re- 
cent progress in crystal growth techniques makes it pos- 
sible to fabricate and control several classes of oxide het- 
ero junctions. Since Ohotomo et at. found interface- 
specific conducting states between two different insu- 
lating perovskite transition-metal oxides^i^, one of the 
main interests has been directed toward the electronic 
properties near the hetero-interface, drastically varying 
with respect to the bulk components. One of the best- 
known examples is the LaAlOs/SrTiOa heterostructure. 
Although both LaAlOs and SrTiOs are conventional 
band-insulators, recent experiments have revealed two- 
dimensional superconductivity^, ferromagnetic correla- 
tion^ and even the coexistences of the two states at n- 
type polar interfaces''^. These phases exhibit a distinct 
dependence on film-thickness and growth-condition^ 
and an intriguing electric-field response^i^. Furthermore, 
high carrier-mobility at the interface^iiSrii^ and two- 
dimensional superconductivit y ^^1^^ have been reported 
for LaTiOs (LT0)/SrTi03 (STO) heterostructures. 

Motivated by these experimental findings, many theo- 
retical approaches have been proposed to describe the na- 
ture of quasi two-dimensional electron systems localized 
around the interface^'^^. Particularly the LTO/STO in- 
terface has been extensively investigated due to the chem- 
ical similarity and small mismatch of lattice constants be- 
tween perovskite compounds LTO and STO^^. Okamoto 
and Millis analyzed a generalized Hubbard model in the 
framework of Hartree-Fock approximatio n^ ^ i ^ ^ and dy- 
namical mean field theory (DMFT)2^i2i. The origin of 
metallic behavior at the LTO/STO interface has been 
clarified, and, moreover, possible ferromagnetic and an- 



tiferromagnetic phases have been investigated, also in re- 
lation to orbital ordering^l and their temperature de- 
pendence^^. First principles calculations have been also 
performed to explore the role of the lattice reconstruc- 
tion, and demonstrated several spin and charge ordered 
phases^^'^^ as well as the effects of lattice relaxationiS, 
which modify the interface electronic structure. More- 
over, Riiegg et al examined the electronic and thermo- 
electric transport properties, particularly focusing on the 
characters of quasi-particles in the inhomogeneous lay- 
ered system^'^. 

In this paper, we present a theoretical study of elec- 
tronic and magnetic properties in correlated heterostruc- 
tures of the Mott-insulator (MI) embedded between two 
band insulators (BI), such as STO/LTO/STO struc- 
tures. A sketch of our model system is shown in Fig. [TJ 
According to the previous optical and photoemission 
studies, LTO/STO hetero junction develops a quasi- two- 
dimensional electron gas extended only a few unit-cells 
across the interface. As mentioned in Ref. [28|, however, 
some experimental setups may allow us to tune the di- 
mensionality of the interface electron system. There- 
fore, reflecting the broadness of the conducting elec- 
tronic state, the possible d band occupation would vary 
abruptly or gradually from n = 1 {d^) at LTO to n = 
{d^) at STO region beyond the interface if the thickness 
of the heterostructure along the stacking direction is large 
enough to reproduce bulk properties in each region. Such 
sharpness of the spatially changing electron density might 
even modify the electronic structure only near the inter- 
face. 

Here we shed light on the nature of the electronic state 
around the hetero-interface with particular emphasis on 
the interplay of charge distribution and magnetic prop- 
erties. As mentioned above, we consider a BI/MI/BI 
heterostructure, taking the example of a LTO/STO het- 
erostructure. For simpHcity, we assume perfect matching 
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FIG. 1: Schematic view of the BI/MI/BI-type [001] het- 
erostructure defined by N layers of positively charged cations. 
We choose the z-axis along [001] direction. Black sites at 
Ri = a{li,mi,ni) correspond to the B sites of perovskite 
ABO3 and A' BO3 lattices, where the relevant conduction 
electrons reside. The A and A' ions located at R\^^ = 
a{li + 1/2, rrii + 1/2, rii + 1/2) are replaced by +1 and neutral 
counter-ions, respectively. 



of the two lattices and, thus, only focus on the electronic 
reconstruction based on the mobile electrons in the Ti- 
derived band. The corresponding model heterostructure 
is described by the generalized Hubbard model defined in 
Ref. [21I. Here our focus is put on the competition of spin 
and charge degrees of freedom while neglecting the or- 
bital degeneracy of Ti-3d band. We treat on-site interac- 
tion with the Hartree-Fock approximation, and calculate 
the ground state phase diagram as a function of the local 
and long-range Coulomb interaction. We obtain a vari- 
ety of the interface-specific magnetic and charge-ordered 
phases. We find that the intimate spin-charge coupling 
in magnetic fields causes a first-order metamagnetic tran- 
sition as well as a reentrant charge-ordering transition. 

This paper is organized as follows. In Sec. [Ill we in- 
troduce the model and method employed in this study. 
In Sec, mil we discuss the T = phase diagram obtained 
as a function of local and long-ranged Coulomb interac- 
tions, and provide detailed studies of some ordered phases 
specific to the interface. Sec. IIVI presents the effects of 
external magnetic fields, and clarifies the origin of field- 
induced spin and charge transitions. Finally, a brief sum- 
mary is given in Sec. |Vl 



II. MODEL 

We study a BI/MI/BI-type heterostructure based on 
the [001] stacking of Mott insulator ABO3 and band 
insulator A^BOs with cubic perovskite lattice structures; 
hereafter we choose the z axis along the [001] direction. 
Taking an example of LTO/STO multilayer structure, A 
(A^) ions represent a valence of +2 (+3) charged cation, 
and mobile conduction electrons reside on the B ions 
which are common constituents for ABOs and A'BOs. 
A schematic view of this system is illustrated in Fig. [H 
We assume the chemical similarity and perfect lattice 



match between two cubic perovskite components. There- 
fore, simulating the different valences between A and A' 
ions, we define the model heterostructure via N counter- 
ion layers with a positive charge +e sitting in the Mott- 
insulating material; A^~^ and A'^^ cations are here substi- 
tuted with +1 and neutral point charges, respectively^!. 
Hence, in addition to the repulsive electron-electron in- 
teraction, electrons suffer from the electrostatic potential 
of the positively charged ions placed between electrically 
active B sites. The charge neutral condition requires 
that the total number of these electrons equals to that of 
+l-charged A site ions. In the present work, the corre- 
sponding microscopic model is simplified with disregard- 
ing orbital degeneracy of B ions, which are sometimes as- 
sociated with interesting phenomena in transition-metal 
oxides. Thus, this model heterostructure captures the 
essential aspects of the competition between spin and 
charge degrees of freedom in a spatially non-uniform en- 
vironment. 

To describe electron dynamics in the background of 
the non-uniform attractive potential, we introduce the 
generalized Hubbard model given by, 

H = Hijand + Hee + ^ei- (1) 

Here, the kinetic term H^and is presented by the single 
band tight-binding model, 

Hband = i4a^jcT + h.c), (2) 

<i,j>,cr 

where Qct is an annihilation operator of an electron with 
spin a =t,i at the site labeled by index i diS Ri = 
a{li,mi^ni) with lattice constant a. Electron hopping 
with the transfer integral t is limited between the nearest- 
neighbor sites on the square lattice. H^e and H^i corre- 
spond to electron-electron and electron-ion Coulomb in- 
teraction given by, 

Fee = UynM + - y (3) 

2 ^ 

H — ^ f4') 

where nia = 4cr^*cr is the spin dependent occupation 
number, and Rf"^ = a{li + 1/2, + 1/2, + 1/2) spec- 
ifies a position of the cation. Both the electron-electron 
repulsion and the electron-ion attraction are introduced 
as a long-ranged Coulomb interaction whose characteris- 
tic energy scale is given by the coefficient Ec = e^/ae 
with the dielectric constant e of the host lattice. In 
Eq. (j3j), the strong on-site repulsion between electrons 
is considered via an ordinary Hubbard interaction term 
with magnitude U. Note that the strength of the on-site 
interaction has the same value on all B sites through- 
out the heterostructure, i. e. the spatial inhomogeneity 
of this model is characterized by the H^i term. There- 
fore, it is naturally assumed that the electron density 



3 



profile is mainly determined by these long-range interac- 
tions. On the one hand, the broad charge distribution, 
extended into the band-insulating material, is favored in 
order to suppress repulsive electron-electron interaction 
i^ee, but on the other hand, the attractive electrostatic 
potential Hei confines the conduction electrons in the 
Mott-insulating region. 

We also take the effects of external magnetic fields into 
account to give a systematical study of the magnetism in 
the inhomogeneous system. In the current study, the 
magnetic fields are treated through the Zeeman term. 

Hz = - ]^9I^bH ^ rriiz , (5) 

i 

where g and jiB indicate the electron g-factor and the 
Bohr magneton, respectively, and rriiz = rii^ — rii^ is a 
magnetization along the applied magnetic fields. Here- 
after, the applied magnetic fields are measured in units 
of Ig/J^B- 

The ground state phases are governed essentially by 
three parameters: thickness of positive ion layers on- 
site Hubbard U and the coupHng constant of long-range 
Coulomb interaction Ec. As mentioned in Ref.[2l|, for the 
example of the experimentally fabricated LTO/STO het- 
erostructure substantial uncertainties exist in estimating 
the realistic values of U and Ec except for an experi- 
mentally tunable parameter N. Particularly the dielec- 
tric constant e = e'^/aEc in bulk STO samples strongly 
depends on temperatur o^^i^^ and electric fields^^, and 
can even change its value by more than one order of 
magnitude as a function of these quantities. Therefore, 
in addition to discussions on the role of U and Ec 
may be important to gain physical insights into the elec- 
tronic and magnetic properties of realistic heterostruc- 
tures, particularly into the interplay of spin and charge 
degrees of freedom. However, a systematic and exact 
computational study of the model Hamiltonian [Eq. ([1])] 
over a wide range of parameters is practically excluded. 
Hence, in the current study, the on-site repulsion and 
long-range Coulomb interaction are treated by employ- 
ing unrestricted Hartree-Fock and Hartree approxima- 
tion, respectively. Here, for the derivation of the ef- 
fective single-particle Hamiltonian, the on-site Hubbard 
term U Ui^rii^ in Eq. (j3j) decoupled as follows, 

+ (4^c4)(4^Qt) (6) 

where the site-dependent expectation values (...) are cal- 
culated self-consistently. However, not only the energet- 
ically most favorable state, several self-consistent solu- 
tions are typically allowed in mean- field treatments. We 
thus consider a variety of ordered states allowed up to 
two sublattices and select the one with the lowest energy 
after many iteration steps of O(IO^). As compared to 
other many-body methods such as dynamical mean field 



theory, the Hartree-Fock treatment treats electron corre- 
lation effects only with limitations. We, however, believe 
that the approximation captures the qualitatively cor- 
rect tendency towards certain electron correlations, while 
overestimating the realization of actual order. 

In the following, we restrict our investigation to the 
physically relevant parameter regime U > Ec zero 
temperature. We consider a supercell geometry of about 
25-layers structure including N cation planes, and em- 
ploy open boundary conditions for solving the mean-field 
Hamiltonian. Since electrons are well confined within 
this multi-layer structure for all the parameter values 
considered in our study, the boundary condition would 
not change the quantitative aspects of our present calcu- 
lation. 



III. PHASE DIAGRAM 

In order to gain insights into the electronic and mag- 
netic properties, we first discuss the ground- state phase 
diagram. In this section, we fix the external magnetic 
field to zero and investigate possible phases in the U-Ec 
plane. As pointed out by previous studies^^'^^, the pro- 
file of the charge distribution is sensitive to the energy 
scale of long-range Coulomb interaction Ec- At small 
Ec^ electrons are extended to the band-insulating region 
due to the weak attractive potential of cation layers, and 
the electron density profile gradually varies from n 1 
(MI) to n (BI) across the interface. In contrast, elec- 
trons are strongly confined along the z direction for large 
Ec^ and a more abrupt change of charge distribution is 
observed in the vicinity of the interface. To character- 
ize this "sharpness" in the charge distribution around 
the interface, we introduce a Thomas- Fermi screening 
length which roughly corresponds to \tf ^ yt/EcC^. 
With this value, we can distinguish between "broad" 
and "abrupt" variations of the charge distribution via 
Xtf < CL and Xtf > cl^ for which the positive electro- 
static potential confines mobile electrons between the in- 
sulating regions. We expect that the electronic and mag- 
netic phases near the interface are sensitive to the elec- 
tron density profile, and change their character around 
Xtf — ci. Similar insights have been also given in the 
preceding studies using double exchange— and Holstein- 
Hubbard modeP^. 

We investigate possible in-plane magnetic and charge 
ordered states with a two-sublattice structure. DMFT 
calculations performed by Okamoto and Millis^^ showed 
that a ferromagnetic ordering is favored in the vicinity 
of the interface in the large U regime. Therefore, besides 
the conventional paramagnetic (PM), ferromagnetic (F) 
and antiferromagnetic (AF) states, we allow for a spatial 
modulation of the magnetic moment. Along the z-axis, 
a possible magnetic state can vary smoothly from the 
AF state in the Mott-insulating region to the FM state 
around the interface. We search for the favorable direc- 
tion of the FM moment relative to that of the AF mo- 
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FIG. 2: Ground- state phase diagram as a function of the on- 
site repulsion U and long-range Coulomb interaction Ec. The 
number of the counter-ion layers is = 4. Our focus is put 
on the physically relevant regime, U > Ec. Six phases are 
shown: paramagnetic metallic state (PMM), paramagnetic 
charge ordered state (PMCO), antiferromagnetic state with 
(7r,7r,7r) spin configuration (GAF), (0, 0, 7r)-type antiferro- 
magnetic state (AAF), GAF state with canted antiferromag- 
netic state at the interface (CA-GAF) and GAF state with 
checkerboard charge-ordering at the interface (CO-GAF). The 
solid and broken lines indicate the first- and second-order 
transitions, respectively. 



ment. In addition, although a simple cubic lattice with 
band filling 1 prefers the G-type AF (GAF) ordering 
showing (7r,7r,7r) spin configuration particularly in the 
N ^ oo limit, the strong spatial inhomogeneity in the 
charge density profile may induce a uniaxial distortion 
for the possible spin alignment in the z-direction. Hence, 
as trial magnetic states in the Mott-insulating region, we 
further take into account the A-type AF (AAF) and C- 
type AF (CAF) structure having the (0, 0, tt) and (tt, tt, 0) 
spin orderings, respectively. 

Fig. [2] presents the calculated phase diagram for an 
A = 4 layer structure with counter- ions at z = ±1.5 and 
z = ±0.5. To obtain the phase diagram, we compare 
the energies of several different ordered phases and se- 
lect the one with the lowest energy. Six phases are found 
within the current mean- field analysis: PM metallic state 
(PMM), PM charge ordered (CO) state with checker- 
board pattern (PMCO), GAF state, AAF state, canted 
AF state around the interface accompanied with GAF 
ordering in Mott-insulating region (CA-GAF), and fer- 
romagnetic checkerboard CO state at the interface plane 
with GAF ordering for Mott-insulating side (CO-GAF). 
The phase boundaries, represented by the solid (broken) 
lines, correspond to the first- (second-) order phase tran- 
sitions. In order to examine the role of the counter-ion 
layer thickness A, the thicker heterostructure, involving 
A = 10 positive layers, is also investigated. The cor- 
responding phase diagram is, however, analogous to the 
A = 4 case apart from the AAF state, which is absent 



in the case of A = 10. This type of A-dependence is 
further discussed in Sec. IIII A[ 

Intriguingly, the charge (magnetic) properties in the 
mean-field phase diagram are dramatically changed 
around Ec/t 1 (U/t c:^ 3). Let us first focus on the 
small Ec regime (weak confinement regime \tf > cl)- At 
small U the possible ground state is PMM. With increas- 
ing f/, this state undergoes a continuous phase transition 
to the GAF state, whose behavior is in good agreement 
with Ref. 113; in spite of the lower occupancy n ~ 0.4, 
the finite AF moment, aligned in the same direction as 
that of the neighboring layers, is observed at the interface 
layer (see also Fig. [3|). Note that this GAF state is sta- 
bilized in the whole Ec region at intermediate U . When 
U is further increased, we encounter a distinct magnetic 
phase: a ferromagnetic component of the magnetization 
spontaneously arises around the hetero-interface, giving 
rise to a tilt of the AF moment. We, thus, arrive at 
the CA-GAF phase, although AAF states are still sta- 
ble under certain conditions. The detailed study for the 
CA-GAF state is given in Sec. IIII Al 

On the other hand, the charge sector exhibits quite 
different properties for Ec/t > 1 {Xtf < cl)- In this case, 
the PMCO state is stabilized in the small-[/ region, which 
turns into the GAF state around U/t 3 via a first-order 
transition. Here, the charge-ordering is realized in the 
whole Mott-insulating region. We naturally expect that 
the PMCO state is stabilized by the long-range Coulomb 
interaction^. For the large U region, we propose another 
interesting mechanism to stabilize the charge ordering 
in the CO-GAF phase, which will be addressed later in 
Sec. UHBl 

In the following subsections, we give a more detailed 
discussions on the phase diagram with particular empha- 
sis on the two remarkable phases realized in the large U 
region: the canted antiferromagnetic state and the charge 
ordered state. 



A. Canted Antiferromagnetic state for Ec/t < 1 

We start with discussions on the possible ordered states 
for Ec/t < 1 {Xtf > cl)' In this case, the weak at- 
tractive potential makes the conduction electron deeply 
penetrate into the band-insulating region, and, thus, the 
charge distribution near the interface is relatively sensi- 
tive to the values of Ec. Hence, the electronic phases for 
Ec/t < 1 would be connected with the spatial variation 
of the charge distribution through the band-filling de- 
pendence of the typical bulk phases^^^^^. In this regime, 
we find three phases which are common to A = 4 and 
A = 10: PMM, GAF and CA-GAF states. Among them, 
the most salient phase we found is the CA-GAF state, be- 
cause there is no analog in the bulk phase diagram^i^i. 
This state is defined by the z-axis modulation of the an- 
tiferromagnetic spin alignment; around the interface, the 
alternating up and down spins are both tilted in the same 
direction, and, thus, a ferromagnetic moment perpendic- 
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FIG. 3: Layer- dependent electron density (filled square), mag- 
netization (open circle) and staggered magnetization (filled 
circle) for the parameters U /t — 12, Ec/t = and = 10. 
Counter-ions are placed z — ±0.5, d=1.5, • • • zb 4.5. The in- 
terface layers are located at z — ±5, separating Mott- and 
band- insulator regions. 



ular to the AF moment emerges. To obtain the above 
solution we consider the two kinds of the phases accom- 
panied with finite magnetization around the interface, 
which are distinguished by the angle between the direc- 
tion of the ferromagnetic and antiferromagnetic moment: 
perpendicular or parallel. With the energetic compari- 
son, we find the solution with perpendicular one is fa- 
vorable, and therefore the interplay of the F and AF or- 
derings induces the spin canting in the vicinity of the 
interface. 

In Fig. [3l several quantities in the CA-GAF state are 
plotted as a function of the distance z from the center of 
the heterostructure: the spatial modulation of the layer- 
dependent electron density n{z)^ magnetization m{z) and 
staggered magnetization nistagiz)- This figure confirms 
the above statement that the weak charge binding Ec < t 
{Xtf > o) induces the charge leakage toward the band- 
insulating region; the electron density drops from n c:^ 1 
to :^ over a few layers across the interface, and the 
finite electron density is even observed in outer layers 
(|2;| > 5). In addition, the band occupation is always 
less than n = 1 through the heterostructure. These 
features are analogous to the preceding studie^iiii^i^^i^. 
Turning to the properties of magnetism, the values of 
the staggered magnetization are almost saturated inside 
the MI region (\z\ < 5) and at the interface layers, re- 
spectively, and the magnetic moments are gradually di- 
minished away from the corresponding regions. We note 
that in this CA-GAF phase, the characteristic behaviors 
of n(z), nistagiz) and m{z) presented in Fig. [3] are con- 
sistent with the DMFT results^^, where good agreement 
between DMFT and Hartree -Fock analysis is only found 
for the electron density and the staggered magnetization. 

All such features manifest themselves in the corre- 
sponding layer-resolved single-particle spectral function 
A(j^cr'{z^uj) in Fig. SI where we plot one of the spin- 
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FIG. 4: Layer-resolved single-particle spectral function A^a 
for the CA-GAF state with U/t = 12, Ec/t = 0.8, N = 10. 
Solid and broken lines indicate A^^ and A^^ components, re- 
spectively. 



dependent components Aa^a{z^uj). Around the center of 
the heterostructure {z 0), the conduction band gets 
quite narrow, and the opening of the gap (c^ U mstagi^) — 
U) is observed due to the AF ordering; inside the MI re- 
gion, the competition of the on-site U and the attractive 
Coulomb potential sustains the electron occupancy ap- 
proximately 1 per site. In contrast, in the BI regions, 
the form of spectral function well reproduces that of 
the free electrons in the bulk system a^^i^^i^^ , and the 
whole conduction band is finally pushed above the Fermi 
level for \z\ > 7. However, when approaching the in- 
terface, the spectral function is gradually shifted down- 
ward to the chemical potential because of the interplay 
of the electron spreading and binding effects along the 
2:-axis. Metallic layers are, thus, formed around the in- 
terface, and as can be seen in Fig. [3l these layers also 
carry a ferromagnetic spin polarization. Note that al- 
though the layer-dependent magnetization looks in good 
agreement with the previous DMFT analysis^^, the ori- 
gin of the present magnetization should be understood 
by the so-called Stoner mechanism. Actually the mean- 
field phase diagram of two-dimensional Hubbard model 
well describes the emergence of the spontaneous mag- 
netization for U/t 10 at the interface showing 0.4 
band occupation?^. The small, but finite ferromagnetic 
moments in the neighboring layers could, thus, be iden- 
tified as a proximity effect from the F ordering in the 
interface. Therefore, it is necessary to employ other so- 
phisticated methods to give more accurate descriptions of 
the magnetism in the correlated heterostructure. These 
treatments are, however, beyond the scope of this paper, 
and we leave this issue for future research. 

We now return to the discussion on the CA-GAF state. 
As can be seen in Figs. [3] and HI the gradual suppres- 
sion of the spin canting is observed toward the center of 
the heterostructure, and accordingly the canted AF state 
is restricted to only a narrow region around the inter- 
face. Therefore, these results suggest the weak thickness- 
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dependence for the stabilization of the CA-GAF state. 
This finding is supported by the obtained phase diagram 
with on-site U and cation layer thickness which is 
presented in Fig. [5l For N > 2, the general properties 
are similar; there are PMM, GAF and CA-GAF states 
from small to large U regime. In contrast, the result 
of = 1 shows a considerably different behavior: the 
PMM state for small /7, the AAF state for narrow range 
of intermediate U and the F state for large U. Here, 
in the CA-GAF state at A/" = 2,3, all layers in the MI 
show the metallic density of states and a canted spin 
state, while such layers are limited around the interface 
for thicker heterostructures A" > 4. Compared with the 
PMM-to-GAF transition line, the phase boundary be- 
tween the GAF and the CA-GAF state are less sensitive 
to A^. As mentioned above, this tendency may be ex- 
plained by the Stoner concept for metallic ferromagnets, 
because the electron occupancy at the interface layer is 
essentially independent of the layer thickness for large 

In contrast, if Ec/t becomes smaller it is expected that 
the corresponding phase diagram shows a behavior more 
sensitive to the thickness N. This point is understood as 
follows; the weak charge binding reduces the occupancy 
at the center of heterostructure, and accordingly smaller 
N further assists the electron leakage toward the band- 
insulating region and shows smaller occupation even at 
the central layers, as compared to thicker heterostruc- 
tures. The calculated U-Ec phase diagram in FigE] con- 
firms the above statement. Actually, at Ec/t 0.4, we 
can find the AAF phase at A' = 4, which is not found 
at N = 10. Unfortunately, however, it seems difficult 
for the present mean-field analysis to catch the physi- 
cal aspect for such small Ec region because the Stoner 
mean field treatment likely overestimates the ferromag- 
netic tendency in those layer s^^a^i^. Hence, the study 
for such Ec/t <^1 regime will be given elsewhere. 



B. Charge ordered state for Ec/t > 1 

In this subsection, we focus on the regime Ec/t > 1 
{Xtf < o) where electrons are strongly confined along 
the z direction. Compared to Ec/t < 1, the electronic 
phase diagram for Ec/t > 1 is less sensitive to the layer 
thickness N and the profile of the electron density. In the 
corresponding parameter regime, the strong attractive 
potential traps the electrons inside the MI region, and 
accordingly the most mobile electrons are absent in the 
band-insulating region. This observation suggests that 
the role of the interlayer hopping, i. e. electron hopping 
along the z direction, is distinguished from that of in- 
planes, particularly around the interface. It is because 
the layer-dependent electron occupancy shows an abrupt 
change across the interface: from n 1 in the Mott- 
insulating region to n in the band-insulating region. 
In the region of Ec/t > 1, the corresponding phase di- 
agram in Fig. [2] typically consists of the COPM, GAF 
and CO-GAF state. Note that the stabilization of the 
COPM phase as weh as GAF state is well understood 
via the bulk phase diagram^^. So, we hereafter focus on 
the interface-specific CO-GAF state. 

The sublattice-dependent charge distribution for the 
CO-GAF state is shown in Fig. [6l Most of the elec- 
trons reside inside the heterostructure, and, thus, only 
a few electrons are found in the outer layers. In ad- 
dition, the resulting checkerboard-type charge ordering 
can be found only at the interface layer, which is accom- 
panied by a ferromagnetic spin configuration. Therefore, 
the present result is consistent with the previous density 
functional calculation performed by Pentcheva and Pick- 
ett^^. The corresponding layer-resolved spectral function 
A{uj) is presented in Fig. [71 where A{A) (solid lines) and 
74(B) (dotted lines) indicate the resulting spectral func- 
tions for the A and B sublattices, respectively. Due to 
the interface charge-ordering, a gap is formed at the in- 
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FIG. 7: Sublattice-dependent layer- resolved spectral function 
of CO-GAF state at U/t = 11, Ec/t = 1.7 for the = 4 
heterostructure. Solid and dashed line represent the spectral 
functions for A- and B-sublattices, respectively, and a i 
illustrate the corresponding spin state for each sublattice. 



terface layer = 2), which would suppress the metallic 
interface state with approaching Ec/t 2.4, where the 
insulating interface is obtained as in Ref. [26|. 

We now examine possible driving forces for the inter- 
face charge-ordering. Figure [8l^a) shows the order pa- 
rameter \nA — tibI of the ferromagnetic CO phase as a 
function of Ec. Here, ua {^b) represents the interface 
electron density at the A (B) sublattice. The filled and 
open circles show the results for different values of in- 
terlayer hopping tz- At t^/t = 0, the electron motion 
is restricted in the two-dimensional sheet while tz/t = 1 
corresponds to the isotropic one. We find that the in- 
terface charge-ordering is induced in both cases. Note 
that the charge-ordering is stabilized even at t^/t = in 
the large Ec region, implying that the ferromagnetic CO 
state is stabilized due to the long-range Coulomb repul- 
sion within the same planed Namely, the interface layer 
has the electron density n ^ 0.4, which increases with 
increasing Ec and finally approaches the quarter filling 
where the checkerboard charge-ordering is particularly 
favorable. 

On the other hand, we should take into account the 
effects of electron transfer in the z-direction to explain 
the charge order emerging for tz/t = 1 in the region of 
Ec/t < 1.5 in Fig.[8fa). Figure [8fb) presents the detailed 
study of the t^-dependent CO parameter at Ec/t = 1.4. 
As seen from the figure, the CO phase becomes more pro- 
nounced with increasing tz- This trend seems, however, 
to contradict the above-mentioned mechanism, since the 
increase in tz causes the electron transfer into the band- 
insulating region and, thus, depletes the interface, po- 
tentially suppressing charge ordering. Therefore, it seems 
difficult to attribute the origin of the CO-GAF state only 
to the long-range Coulomb repulsion for Ec/t 1.5. 

To elucidate another possible mechanism, we focus on 
a nontrivial role of tz that induces the effective inter- 
layer magnetic interaction. Figures M,^) and (b) show 
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FIG. 8: Order parameter \nA — nB\ of the interface charge- 
ordering for the N — \ heteroj unction with interface electron 
occupancy at A(B)-sublattice nAins)- (a) \nA — nsl is il- 
lustrated as a function of Ec for U/t — 11. Open and filled 
circle denote the resulting isotropic {tz/t = 1) and anisotropic 
(tz/t — 0) electron transfer, respectively, (b) the detailed tz- 
dependence is depicted for Ec/t — 1.5. To obtain the values 
for tz/t ^ 1 we here use the result of t^/t = 1 as an initial 
input and tune the magnitude of the inter-layer hopping. 



the sublattice-dependent magnetization and the inter- 
face CO parameter as a function of [/, respectively. In 
the vicinity of the CO-transition point U/t c:^ 9.5, the 
interface magnetization spontaneously appears via the 
exchange mechanism, causing different spin configura- 
tions between two sublattices; A-sites (lower occupancy) 
present a ferromagnetic spin alignment with the neigh- 
bor layer located inside the heterostructure while B-sites 
(higher occupancy) prefer the anti-parallel one. Hence, 
ferrimagnetism is stabilized at the interface. According 
to our numerical results, A-sites always carry the lower 
electron density compared to B-sites around z '2:^ 2. To 
clarify the relation between this spin structure and the in- 
terface charge-ordering, we take into account the charge 
confinement effect along the z-axis; for Ec/t > 1 (i.e. 
Xtf < o)^ the strong charge binding via the positive 
background localizes the electrons inside the heterostruc- 
ture {z < |3| in Fig. [6]), and thus it is not energetically 
favorable to keep the higher occupancy outside the het- 
erostructure. Therefore, the electron occupancy for the 
A-sites decreases to suppress the charge-transfer toward 
the band-insulating region. This in turn increases the 
occupancy at B-sites, which may enhance the antiferro- 
magnetic correlation along the z direction. This mech- 
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FIG. 9: Plot of [/-dependence of several quantities near the 
interface around the transition point from GAF to CO-GAF 
state, (a) the magnetization for A-sublattice (squares) and 
B-sublattice (circles) at Ec/t — 1.4. Open and filled symbols 
refer to the values of the interface (z — 2) and the neighbor 
layer {z — 1), respectively, (b) corresponding order parameter 
of CO state at interface layer. 



anism based on the interlayer spin couplings naturally 
explains the behavior of the CO parameter for small tz 
in Fig.[8](b). In summary, we may interpret the interface 
charge order as being induced by interlayer spin coupling 
in order to reduce the charge leakage from the interior 
toward the band insulator. This spin-driven mechanism, 
together with long-range Coulomb interaction, stabilizes 
the CO-GAF phase even for the relatively weak Ec. 



IV. EFFECTS OF MAGNETIC FIELDS 

In this section, we discuss the nature of the het- 
erostructure under external magnetic fields to further 
clarify the interplay of the magnetism and the charge 
density profile. As discussed in the previous section, the 
possible ground states are closely connected with the pro- 
file of charge distribution. At zero magnetic field, the 
long-range Coulomb interaction Ec plays an important 
role in determining the charge density profile along the 
z-axis. We also note that magnetic fields can have a 
strong influence on the charge distribution by modifying 
in the spin structure which governs the charge distribu- 
tion in the vicinity of the interface^^. In the following, 
we analyze the effect of magnetic fields on the density 
modulation. Throughout this section, we consider the 
heterostructure with = 10, but the obtained results 



apply to a wider range of N (see Sec. IIV Ap . We focus 
on the relatively large U region in which the GAF state 
is stabilized inside the heterostructure, and consider the 
magnetic fields perpendicular to the AF-moments using 
the Zeeman coupling term. Under these circumstances, 
we will deal with three phases, i. e. GAF, CA-GAF and 
CO-GAF. 



A. Metamagnetic transition 

We first concentrate on the appearance of a meta- 
magnetic transition in the heterostructure under external 
magnetic fields. Figure [TOT a) shows the layer-dependent 
magnetization as a function of H for fixed values U/t = S 
and Ec/t = 0.8. Note that in the absence of external 
fields, the corresponding ground state is the GAF in this 
parameter regime (Sec. lIIip . With increasing H the inter- 
face magnetization (open circles) displays a discontinuous 
jump around H ~ 0.8, signaling a first-order metamag- 
netic transition. This is confirmed by a hysteresis loop 
emerging with rising and lowering magnetic fields. In 
contrast, we do not find any anomalous behaviors in the 
magnetization at the central layer (filled circles); it in- 
creases monotonically without exhibiting any jumps, and 
smoothly approaches the saturated value m c:^ n 1. 
These results imply that the metamagnetism observed 
here is intrinsic to the interface. Interestingly, this meta- 
magnetic transition is intimately connected with the 
charge distribution, as shown in Fig. [TOT b). where the 
left (open) and right (filled symbols) axes display the 
computed electron density at the interface and the cen- 
tral layer, respectively. As clearly seen in this figure, 
the electron density exhibits the same anomalous behav- 
ior as in the magnetization in Fig. [TOT a): (i) there exist 
the abrupt jumps in both quantities only at the inter- 
face at the same magnetic fields and (ii) a hysteresis loop 
emerges with increasing/decreasing magnetic fields. This 
pronounced change in the electron density is restricted to 
the vicinity of the interface. 

Therefore, in order to understand the origin of the in- 
terface metamagnetism, it is crucial to elucidate how the 
magnetism is related to the charge distribution in the 
heterostructure. To this end, we show the [/-dependent 
electron density for H = at several layers in Fig. [11) 
the center layer {z = 0), the nearest-neighbor layer to 
the interface {z = 4) and the interface layer (z = 5). The 
characteristic feature appears at the PMM-to-GAF tran- 
sition point U/t 3, around which the gradients of the 
curves for z = 4 and z = 5 change their signs, in con- 
trast to the monotonic change in the PMM state (dotted 
lines). On the other hand, the electron density at the cen- 
tral layer depends only weakly on U. This observation 
is consistent with the previous finding that the density 
profile is affected by the electronic phases inside the het- 
erostructurei^A^. Particularly in Ref. [sH, Riiegg et al. 
suggested that the visible modification occurs via the lo- 
calization of the electronic states with undergoing a Mott 
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FIG. 10: (a) Comparison of magnetization curves at the in- 
terface (open) and at the center of the heterostructure z — 
(filled symbol) with the fixed values U/t — 8 and Ec/t = 0.8. 
(b) The corresponding electron density at the interface (left 
axis) and z = (right axis). The up (down) arrow, denot- 
ing the hysteresis loop, presents the magnetization under the 
increase (decrease) of magnetic fields. All quantities are com- 
puted under the magnetic fields perpendicular to the AF mo- 
ment in the Mott-insulating region. 



transition; mobile electrons are tightly bounded inside 
the heterostructure to keep the higher band-occupation 
c± 1. This scenario is also applicable to the present case 
since in the bulk phase diagram, the conventional GAF 
phase prefers the singly occupied sites. Therefore, in or- 
der to support the emergence of antiferromagnetism, it is 
necessary to confine the electrons inside the heterostruc- 
ture with varying the density profile along the z-axis. We 
note that such charge redistribution occurs only near the 
interface^^, and plays a rather minor role in comparison 
with the long-range interaction Ec. Hence, the interface 
electrons are mainly confined to the Mott-insulating re- 
gion, and the resulting density profile changes abruptly 
across the interface. 

Recall here that the applied magnetic fields have a ten- 
dency to destroy the AF ordering. Such a suppression 
seems to be stronger near the interface than the central 
layers because the number of the nearest-neighbor sites 
coupled antiferromagnetically to each other decreases 
from 6 at the central layers to 5 at the interface. This, 
in turn, increases the field-induced magnetic moment 
around the interface. Especially the interface magnetiza- 
tion is immediately saturated because of the extremely 
weak AF ordering and the large spin susceptibility via 
in-plane instability toward Stoner ferromagnetism. Nev- 




U/t 

FIG. 11: The [/-dependent electron density for given layers 
labeled by z = 0,4, 5 at Ec/t = 0.8 and zero magnetic field 
(see Fig. [3|): broken lines for the paramagnetic state at the 
interface (z = 5, lower panel) and for the neighboring plane 
to the interface {z = 4, upper panel). For layer indices, see 
FigEl 

ertheless, below the metamagnetic transition point, the 
interface AF-moment remains finite through the proxim- 
ity effect from the neighboring layers (see Fig. [3]). Thus, 
it will be energetically favorable in a magnetic field to 
increase the electron occupation at the interface. On 
the other hand, the charge confinement along the z-axis, 
driven by magnetism, would become more relevant for 
large U because of the strong AF coupling as shown in 
Fig. [TTl This means that the competition between cor- 
relation effects and the influence of the magnetic field 
is crucial for determining the AF order at the interface. 
This gives rise to the metamagnetic transition accompa- 
nied by a redistribution of electron charge near the in- 
terface. At the critical magnetic field, we indeed find the 
discontinuous drop in the staggered magnetization near 
the interface, confirming that the metamagnetism is in- 
deed caused by the charge redistribution of electrons, as 
seen in Fig. [TTl 

This view is further confirmed by our numerical results. 
Figure [12] shows the computed magnetization at the in- 
terface for several values of U. The resulting magneti- 
zation remarkably changes its behavior across a critical 
end point of U/t 7.5. Though the magnetization varies 
continuously below this value, the first-order metamag- 
netic transition is obtained for [//t > 8, above which the 
critical magnetic field displays a gradual decrease with in- 
creasing U. These tendencies agree well with the above 
suggestion, because relatively large values of U may en- 
hance the difference between two effects: spreading and 
confinement of electrons along the z direction due to mag- 
netic fields and local Coulomb interaction, respectively. 



B. Field induced charge ordered state 

In the relatively large Ec regime, the charge redistri- 
bution under external magnetic fields reveals another in- 
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FIG. 12: The magnetization curves at the interface for several 
values of the on-site interaction. The corresponding ground 
states ai H — are GAF for U /t — A (filled squares), 8 (open 
circles), and CA-GAF dXU /t = 12 (filled circles), respectively. 



teresting aspect of the correlated heterostructure: field- 
induced charge ordering. Generally, the strong inter-site 
Coulomb interaction (c^ Ec) may break the translation 
invar iance and induce the CO in the bulk systems. As 
mentioned in the previous section, at = 0, there ac- 
tually exists the interface CO state of checkerboard pat- 
tern, i.e. CO-GAF phase^^ Figure [TSTa) displays the be- 
havior of such interface charge ordering in applied mag- 
netic fields. In weak magnetic fields, the CO immediately 
disappears with increasing field. This is caused by the 
emergence of the special magnetic order at the interface 
which partially supports the stabilization of the interface 
charge-ordering as mentioned in Sec. IIIIBi so that the 
disappearance of the magnetic order weakens the inter- 
face charge-ordering. 

Surprisingly, we encounter a reentrant CO state with 
further increasing H] the checkerboard charge-ordering 
emerges again around 0.9 only at the interface. Note 
that there is a discontinuous jump accompanied by a hys- 
teresis loop in the order parameter, implying that the 
transition is of first order. Furthermore the order param- 
eter of the field-induced CO phase is almost unchanged 
once it is induced. This i^-dependence suggests that 
the origin of the interface charge-ordering is the inter- 
site Coulomb repulsion. This view is indeed confirmed 
in Fig. [TST b). Above H 0.9, around which the field- 
induced CO phase occurs, the averaged electron density 
at the interface riave = (^a + nB)/2 shows a plateau at 
^ave =0.5. As is the case for a conventional GAF state, 
the checkerboard CO can be caused by the Fermi-surface 
nesting, particularly at quarter filling. Therefore, the 
in-plane CO state stabilized with riave = 0.5 above the 
critical magnetic field in Fig. [T3l evidences that this CO 
state is driven by the intersite Coulomb repulsion under 
the quarter-filling condition. We note in Fig. [TSTb) that 
below the critical magnetic field {H ~ 0.9), there is an- 
other first-order transition around H ~ 0.85 with a small 
jump in Uave- This is due to the metamagnetic transition 
discussed in Sec. IIV Al (see also Fig. [HTa)). 
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FIG. 13: Plot of interface electron density dXU/t — 11^ Ec/ — 
1.5 as a function of magnetic fields, (a) order parameter of the 
checkerboard charge-ordering. The increase (up arrow) and 
decrease (down arrow) of the magnetic fields yield a hystere- 
sis loop, (b) the corresponding interface occupation number, 
defined by riave = {tia + ub)/'^- Dotted line shows the value 
of the quarter filled band: Uave = 1/4. 



Since there exists a finite charge excitation gap at the 
interface, the field-induced CO state would be experi- 
mentally detected, for example, as a rapid drop in the 
conductivity with increasing H. 



C. Phase diagram under magnetic fields 

It is instructive now to present the phase diagram un- 
der magnetic fields. Figure [HTa) shows the calculated 
phase diagram with the parameter Ec and the mag- 
netic field H. The field-induced charge-ordering (closed 
squares) can be found in the regime of relatively large 
Ec and H. This trend indeed confirms our prediction 
presented in Sec. IIV Bl concerning the origin of the tran- 
sition: the important interplay of the band occupancy 
(^ave = 0.5) and the long-range Coulomb repulsion. In 
contrast, increasing Ec seems less relevant for the meta- 
magnetic transition (filled circles), since this transition 
is accompanied with the increase of electron density at 
the interface, which is unfavorable for large Ec. Interest- 
ingly, such charge redistribution would further enhance 
the tendency toward the field-induced CO transition. In- 
deed these transition lines gradually merge in the large 
Ec region. 

The nature of phase transition from CA-GAF to CO- 
GAF state under the magnetic fields, shown by the dot- 
ted line in Fig. HH^a), remains still open. At intermedi- 
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FIG. 14: The obtained phase diagrams under external mag- 
netic fields: (a) H-Ec phase diagram for U/t = 11. Open 
circles: CO- and CA-GAF phase boundary at zero-magnetic 
field; filled circles: the first-order metamagnetic transition 
from CA-GAF to GAF state with paramagnetic interface 
(PM-GAF),; filled squares: first-order phase boundary of the 
field-induced CO state. As well known, the PM-GAF state 
turns into the paramagnetic phase for large H. (b) H-U 
phase diagram with different values of Ec. Filled circles: the 
first-order metamagnetic transition line for Ec/t = 0.8, across 
which GAF (or CA-GAF) turns into PM-GAF with the in- 
crease in H; the characteristic feature is, however, changed 
to continuous one for U/t ^ 7.5. Filled squares: the phase 
boundary obtained at Ec/t = 1.5, illustrating the discontin- 
uous transition toward the field-induced CO phase with in- 
creasing H. 
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FIG. 15: Schematic H-Ec phase diagram with particularly 
focusing on the boundary between CO-GAF and CA-GAF 
phases. Solid lines: first-order transition, broken line: second- 
order transition. Thus, (a) CO- and CA-GAF states are sep- 
arated by a discontinuous transition, while (b) there exists a 
critical point (CP), above which these phases are connected 
continuously. 



(a) and (b), is correct. 

Figure [Mf b) presents the U -dependence of the meta- 
magnetism and the field-induced charge-ordering, respec- 
tively. As mentioned in Sec. IIV A| the metamagnetic 
transition can appear only in the large U region, and 
thus a critical point may exist around U/t 7.5. Here, 
the overall [/-dependence of the field- induced charge- 
ordering resembles that of the metamagnetic transition; 
the value of the critical magnetic field He increases mono- 
tonically with decreasing U. However, the difference be- 
tween these transitions is found in the lowest critical val- 
ues of He. Although the interface metamagnetism tran- 
sition is of first-order only for U/t > 8, we find the dis- 
continuous charge-ordering transition in a much wider 
range of the on-site interaction. This is because the in- 
tersite Coulomb interaction with the specific occupancy 
of electrons plays a dominant role in the transition to 
the field-induced CO phase, and the on-site Hubbard in- 
teraction merely modifies the electron occupancy at the 
interface with inducing magnetism. 



V. SUMMARY AND CONCLUSION 



ate Ec (— 1.4t), the specific interlayer magnetic coupling 
might drive the CO-GAF phase at = [Sec. HITB] . 
Therefore, if such a spin structure is suppressed with the 
increase in H, the CO-GAF state would undergo a first- 
order transition to the CA-GAF state as is the case for 
H = 0. Note that these two phases are distinguished by 
the direction of the interface F-moment as well as the CO 
parameter, \nA — However, for large Ec ^ 1.6t, the 
interface charge-order even survives in the large H re- 
gion without the spin coupling between layers, and thus 
the difference between CO- and CA-GAF states are only 
given by the order parameter \nA — Accordingly, 
as shown in Fig. [15] (b), the CO- and CA-GAF phases 
might be continuously connected beyond the certain crit- 
ical value of H. In the present study, however, we cannot 
determine which of the two scenarios, depicted in Fig. [15] 



In this paper, we have considered the generalized Hub- 
bard model and presented the mean-field study on the 
ground-state properties of the correlated heterostructure 
composed of a Mott insulator (MI) sandwiched between 
two band insulators (BI). The model structure is elec- 
trostatically characterized by the positively charged ions 
located in the Mott-insulating region, causing a spatial 
modulation of the electron distribution. Such a profile 
of electron density along z direction can be tuned by 
the long-range Coulomb interaction and also by the one- 
site interaction via the magnetic moment induced inside 
the heterostructure. Using a relatively wide range of pa- 
rameters and tuning spin structures via external mag- 
netic fields, we have investigated how the inhomogeneous 
charge distribution leads to the intriguing electronic and 
magnetic phases realized in the model heterostructure. 



12 



At zero magnetic field, our model exhibits rich elec- 
tronic phases in the U-Ec plane, including CA-GAF, 
CO-GAF and some other conventional ordered states. 
The resulting phase diagram changes its charge charac- 
ter around Ec/t 1 {Xtf — o)- For Ec/t < 1, where 
the density profile of electrons gradually changes beyond 
the interface layer, we have found the CA-GAF phase for 
relatively large /7, which consists of the GAF oder in- 
side the MI region and a canted AF phase in the vicinity 
of the interface. Interestingly, even though the spatially 
modulated magnetization is caused by the Stoner mech- 
anism, our results are in good agreement with the earlier 
DMFT study of Okamoto and Millis^"^. Note that the 
CA-GAF phase has been proposed for the first time in 
this paper . On the other hand, for Ec/t > 1, we have 
found the ferromagnetic CO state with checkerboard pat- 
tern at the interface as suggested by previous ab initio 
calculations^. In some parameter regime, however, this 
CO state might not be explained only in terms of the 
long-range Coulomb repulsion. As an additional mech- 
anism supporting the interface charge-ordering, we have 
proposed the importance of the interlayer electron hop- 
ping, which causes the site-dependent interlayer spin cou- 
plings and thereby enhances the instability toward the 
charge-ordering together with the long-range Coulomb 
repulsion. 

It has been demonstrated that a strong coupling be- 
tween spin and charge degrees of freedom is crucial for the 
system in magnetic fields, giving rise to the field-driven 
charge redistribution in the vicinity of the interfaces. One 
of the remarkable consequences of this effect is a first- 
order metamagnetic transition. It has been shown that 
this transition is induced by the competition of two ef- 
fects: the localization of electrons inside the center region 
of the heterostructure via the GAF ordering (/7), and the 
charge transfer toward the interface via the spin polar- 
ization {H). This mechanism is fully supported by the 
computed U -dependence of the metamagnetic transition 
and the specific jump in the interface electron density as 
a function of H. 

Another remarkable phenomenon we found here is the 



field-induced charge-ordering transition. This transition 
at the interface is possible only for relatively large Ec^ 
and is explained in terms of the modification of interface 
electron density through the competition of U and H; 
the intersite Coulomb repulsion together with an emer- 
gent quarter-filled band condition gives the instability 
toward the CO phase. This also demonstrates the im- 
portance of the interplay between the spin and charge 
degrees of freedom in magnetic fields. In order to gain 
further insights into the effects of magnetic fields, we have 
presented the Ec-H and U-H phase diagrams, and found 
that characteristic features of the resulting phase bound- 
aries are consistent with our scenario for the field-driven 
metamagnetic/charge-ordering transitions. 

In this paper, we have performed the mean- field anal- 
ysis for the magnetic properties of the heterostructure, 
which are strongly coupled to the inhomogeneous charge 
distribution. Although this simple treatment provides an 
insight into physical properties of the model, it is partic- 
ularly important to confirm the present results by taking 
into account strong electron correlations with many-body 
techniques such as DMFT^^ and slave-boson method^. 
This issue will be studied elsewhere. 
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